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The matching between Schrodinger Functional renormalization schemes and conventional per- 
turbative schemes is usually done using an intermediate lattice scheme. We propose to do the 
matching directly. This requires the perturbative evaluation of Schrodinger Functional correlation 
functions in the continuum. We use dimensional regularization but due to the lack of translational 
invariance in the Euclidean time direction, we employ a general technique introduced by Liischer. 
In this talk I describe this technique and its application to the one-loop expansion of correlation 
functions used in the definition of the renormalized quark mass in the Schrodinger Functional 
scheme. The divergent parts are identified and the computation of finite parts is in progress. 
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1. Introduction 

The QCD Schrodinger Functional (SF) [p], ||] can be used to define finite volume renormaliza- 
tion schemes which are suitable for non-perturbative renormalization. This means that the running 
of renormalized parameters and operators can be traced non-perturbatively from low to high en- 
ergies through numerical simulations. Applications include the running coupling and quark mass, 
moments of structure functions, the static light axial current and four-quark operators (see [Q| and 
references therein for previous applications). The perturbative, high energy matching with con- 
ventional perturbative schemes such as MS, is usually established by using an intermediate lattice 
renormalization scheme. It is possible, however, to do this matching directly through dimensional 
regularization of the SF. This is the aim of our work. Proceeding this way we expect advantages 
for the renormalization of four-quark operators and useful checks of previous results. 

Conventional dimensional regularization techniques rely on the form of propagators in mo- 
mentum space. Due to Dirichlet boundary conditions in the Euclidean time direction the standard 
techniques cannot be applied to the SF. Therefore, we will employ a method due to Liischer JjJ|, us- 
ing the heat kernels representation of propagators. As a first exercise we study the renormalization 
of the quark mass in the SF renormalized scheme [S]. 

2. SF renormalized quark mass 

We consider QCD with two flavors of mass degenerate quarks. The starting point for the 
definition of the SF renormalized mass is the PCAC relation: 

d ll A a ^(x)=2mP a (x), (2.1) 

where A a Ax) and P a (x) are the isovector axial current and pseudoscalar density respectively. This 
relation holds true, up to contact terms, when inserted in any correlation function. The axial current 
is protected against renormalization by chiral symmetry so it stays finite. Thus if we define the 
renormalized mass nip = Z m m through the PCAC relation, its renormalization will be given once 
an independent renormalization condition is imposed on the pseudoscalar density, i.e. Z m = Zp l 
given Pg (x) = ZpP a (x) . The renormalized pseudoscalar density is defined within the SF framework. 

The SF geometry consists in a space-time manifold with Dirichlet boundary conditions in the 
Euclidean time direction and periodic boundary conditions (up to a phase for fermionic fields 1 ) 
in the spatial directions. To prepare the ground for dimensional regularization we consider a D- 
dimensional manifold with an arbitrary number 3 + d of L-periodic spatial directions. Our nota- 
tional conventions for indices and vectors are the same as in External momenta have only 
physical components and we take the boundary fields to be independent of the extra-dimensional 
coordinates. We impose vanishing gauge fields at the Euclidean time boundaries so that the induced 
background field vanishes. 

To define Zp within the SF framework we consider bilinear source fields at the boundaries 
xo = andxo = T : 

= Jd D - l yd D - l zZ(y) Y5 l -T a m, ^' a = Jd^fd^ztim^'®, (2-2) 
1 This phase plays no role in the following discussion and will be omitted to ease the notation. 
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where £, £, £' and C, are the usual boundary quark fields [0] and T a are Pauli matrices acting in fla- 
vor space. Boundary fields renormalize multiplicatively with a common renormalization constant 
Zr 0]- We now introduce the following correlation functions: 



2(D-1) 



(2.3) 



V w " J 3 
where the powers of L are such that the correlation functions are dimensionless. We write their 
perturbative expansion on the renormalized coupling gR as: 

oo oo 

/p(*o) = I f?(x )g 2 R n , fi = £ (2.4) 

Now Zp is defined equating a ratio of renormalized correlation functions to its tree level value[||]: 

#f-f-\ - , . AO). 



fp(xo) R _ ZpZ 2 r.fp(*o) _ v Mxq) _ f P °\x ) 



V(fi)R 



(2.5) 



where correlations are computed in the chiral limit and xo = T/2. To establish the matching be- 
tween the SF and the MS renormalized mass we must consider the perturbative expansion: 



n=0 



(2.6) 



Here we have used Z m = Z p l . In the following section we will see how to compute fp^ and /j 
using dimensional regularization. 

3. One-loop diagrams 

It is not difficult to obtain the perturbative expansion of our correlation functions. At order g\ 
we find the following diagrams: 



f P 1)a (xo) 



f P 1)a '(xo) 



f P 1)b (xo) 



n 



(l)fl 



(\)a> 



(1)6 





*o x 



x' y ' x 




x! Q ,x% x 



o x' Q y ' t o^o *o r JTq 4' T 



Arrows represent free quark propagators and curved lines represent free gluon propagators. 
The pseudoscalar density insertion is marked with a cross and vertical lines are Euclidean time 
boundaries. Self-energy and vertex diagrams are divergent. They can be compactly written like: 

/ dx' dx'o Tr{/(xo,Xo)I(xo,Xo)}, 



V(f) 



dx dx f Q Ti{f(x ,x'o)V(xo,x Q ,XQ)}. 



(3.1) 



Here £(xq,Xq) and Vp(xo,x? ,Xq) are the amputated self-energy and vertex diagram respectively and 
they contain a sum over momenta. Different self-energy diagrams, a and a', are obtained changing 
the external function /(jc^jc^), which is a string of free propagators and gamma matrices. 
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3.1 Dimensional regularization technique 

To regularize the one-loop diagrams we employ a general dimensional regularization tech- 
nique as introduced by Liischer in [Q]. We begin with a suitable representation for free quark and 
gluon propagators similar to the one given in except that here we work in a time-momentum 
representation: 

5(°)(x ,yo;p) = [-d Y -ip- k r k }[P+G ND (xo,y ;p)+P-G DN (x ,yoM 

Dfl(x ,y ;p) = VvoG^o^p) + 8^ k G DD (xv,y Q ;p). (3.2) 

Here P± = i(l ±70) are Dirac space projectors and G B ° Br (xo,>'o;p) are Green functions of the 
extended Laplacian operator, A = — d$ + p 2 , with Dirichlet (D) or Neumann (N) conditions at 
xo = and xo = T. The advantage of this representation is that now we can express the Green 
functions in terms of heat kernels: 

poo 

G B ° BT (x ,y Q ;P)= dtK? QBT (x ,y ;p), (3.3) 
Jo 

where K?° Bt (xo,yo;f>) satisfy the heat equation: 

(d t + A)K BoBT (xo,y Q ;p)=0. (3.4) 

Contributions from different dimensions factorize, in particular we can take out a factor which 
depends on extra, non-physical, components of the momentum: 

K BoBT (x ,y ;p) =K?° BT (x ,y ;p) ]> ^. (3.5) 

Here K? oBt (x ,yo>p) is the heat kernel of the Laplacian in 4 dimensions with the given b.c. 

For illustration we will consider the vertex diagram. Factorizing the contributions from extra 
components of the momentum it can be written in the form: 



poo 

V(f)= / dt l dt 2 dt 3 [r(t 1 +t 2 +t 3 )} d I(f;h,t 2 ,h) 
Jo 

00 

with: r{t)=LT l £ g-'^W p(n) 



inn (3-6) 



It is at this point that we can take the number of extra dimensions d to be a complex number 
d = —2e, and the diagram becomes a meromorphic function of e. At £ = the integration diverges 
due to singularities in the integrand I{f',h,t2,h) when all proper times go to zero simultaneously. 
Since the function r(t) goes as (47Tf) -1 / 2 for small t, the prefactor \r(t\ +^2 + ^3)] 2e acts as a 
regulator and cancels the divergences of the integrand when Re(e) is large enough. The diagram at 
£ = must be defined by analytical continuation as we will see in the following subsection. 

3.2 Analytical continuation to e = 

We now look in some detail at the integrand: 

Kf\h,h,h) = -c f l 3 [ T dx' Q dx'^e-^ +t ^ 2 Tv{ Yk [r^< D (4i„X(i„,^)- 

JO p 

-p 2 r B K» D (x> 0l x )K» D (x ,4)} K° D (xl4) 7k f(x> ,4)} + ■- (3-7) 
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Here Cp = (N 2 — l)/2N with N the number of colors. Omitted terms have the same structure. 
We have factorized the kernels for the 3 spatial directions and the one-dimensional kernels for the 
Euclidean time direction: Kf° BT (xQ,yo). An explicit expression for them can be found in |^|. F A 
and F B are some Dirac structures. To analytically continue to e = we first must identify at which 
values of proper times this expression is not integrable in the absence of regulator. Singularities 
arise at small values of proper times. We study the asymptotic behavior of the integrand when one, 
two, or all proper times go to zero and find that, as anticipated in the previous section, only in the 
latter case we obtain a non integrable expression. It is convenient then to do the following change 
of variables: t\ = ts\, h = ts%, t 3 = ts 3 , with si + s 2 + s 3 = 1 and the vertex diagram can be written: 
poo ,.y 

V(f)= dt t 2 [r(t)}- 2£ ds l ds 2 ds 3 8(s 1 +s 2 + s 3 -l)I{f,ts 1 ,tS2,ts 3 ). (3.8) 
Jo Jo 

To analytical continue to e = we need to know the asymptotic expansion of I(f,ts\,tS2,ts 3 ) 
including terms of order ?~ 3 . Since exponentially suppressed terms do not enter in the asymptotic 
expansion, we can consider reduced heat kernels (see [^). Furthermore, for the vertex diagram it 
is found that only the free space contribution, i.e. the one we would obtain if the boundaries were 
not present, contributes to the divergence. We can then, from the very beginning of the calculation, 
split propagators and the corresponding kernels in its free space part (f) and a surface part (s) [^] : 

S = S f + S\ D A v=^-+D^. (3.9) 

Contributions to the diagram with any surface propagator are finite and can be computed without 
introducing the heat kernel representation. For the contribution coming from considering the free 
part for all propagators we use heat kernels and find: 

P \f;tsi,ts 2 ,ts 3 ) t ~°A(f;s u s 2 ,s 3 )r\ (3.10) 

A{f;s\ ,S2,s 3 ) is a functional of the external function /, well behaved in the domain of integration. 
The superscript / reminds us that we are working only with free kernels. Since the leading term 
goes like t~ 3 , it is not necessary to go further in the asymptotic expansion. Inserting this expression 



in Eq. (3.8) and using the asymptotic behavior of the regulator, we perform the integration over 
proper times at the lower end of the domain of integration and obtain the divergence as a pole in 
e. To make the analytical continuation to e = we add and subtract this divergence and obtain the 
following representation of the diagram: 



e 



V f {f) £ =° J o ds\ ds 2 ds 3 8(\-s\-S2-s 3 ) ^A(f;s u s 2 ,s 3 ) 

+ J o dtt 2 [lf(f,ts u tS2,t S3 )-e(l-t)A(f;s u s 2 , S3 )t- 3 ]^ + 0(s). (3.11) 

The first term is the divergence and the remaining terms are finite contributions. Proceeding anal- 
ogously with the self-energy diagram we obtain the following divergences: 



(i) _ -3C F (o) 

fl '^""sW 71 • (3 - 13) 
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The first term in the sum comes from the self-energy and the second from the vertex diagram. 



Introducing these results in the expansion of Eq. ( |2.q ) we find the correct divergence for the quark 
mass: 

pole — T"7 T • (3-14) 

4. Conclusions and final remarks 

We have seen how to apply dimensional regularization to SF correlation functions. In particu- 
lar we have identified the divergences of the one-loop expansion of the correlation functions used 
to define the SF renormalized quark mass. However it is the finite part what is required to establish 
the matching with other perturbative schemes. 

The computation of the finite part can be divided in different pieces. First of all, the free gluon 
propagator has a particular expression for zero momentum. We must compute the zero momentum 
contribution to our diagrams separately then. This has been obtained analytically using MAPLE. 
On the other hand, for the computation of the divergence only reduced kernels were required, 
in particular for the vertex diagram only the free part contributes. We can compute the remaining 
contribution working with S s and D%%, which have a simpler analytical expression. Here only a sum 
over momenta and two integrations are required. Finally we have the finite contributions coming 



from Eq. ( |3. 1 1| ). Obviously the difficult part here comes from the subtraction of the divergence 
from the integrand. There are five integrations and the sum over momenta to be done. All these 
computations are work in progress. 
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